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The values of four-quark operators averaged over pions are expressed through those averaged 
over vacuum. The specific values are obtained in the framework of the factorization assumption. 
For the condensates of the light quarks of the same flavour qFqqTq the scalar condensate is shown 
to be an order of magnitude larger than the other ones. The condensates containing the strange 
04 ■ quarks qqss appear to be only about twice smaller than those of the light quarks. The degeneracy 

' of the ground state in the Nambu-Jona-Lasinio model is shown explicitly. 

<N 

£ ' I. INTRODUCTION 

3 ' 
^_ 

The matrix elements of four-quark operators in hadronic states provide information on the structure of the quark- 
antiquark sea in the hadrons describing the correlations of qq pairs. Under certain reasonable assumptions the 
expectation value of the scalar quark operator qq was found in jfj]- || to be equal (up to normalization factors) to the 
■ total number of quarks and antiquarks in the hadron. The known expectation values of qq in nuclcons and pions 
W correspond to njy ~ 8 and fi T w 12, exceeding the number of the valence quarks. This signals the large numbers 
ON of the sea quarks even at low energies. 

Another application is the QCD sum rules in nuclear matter. In-medium values of the four-quark condensates is 
an important ingredient of the approach |^|,^| . In the gas approximation it is expressed through the expectation value 
in the nucleons. 

As it stands now, the only calculation of the four-quark condensate in hadrons is that carried out by Celenza et al. 
. In 0| the condensate in the nuclcons was obtained in the framework of the Nambu and Jona-Lasinio model. 
Some data on its value in nucleons were obtained by Johnson and Kisslinger || by considering the QCD sum rules 
for nucleons and isobars. 

I 1 Here we present calculation of the expectation values of the four-quark operators in pions. Besides providing 
Oh, information on the structure of pions, the analysis is the intermediate step in the calculation of the expectation values 
in nucleons in the framework of models where the quark-antiquark sea is contained in the pion cloud. Some of the 
models analyzed in || as well as perturbative chiral quark model JTg) ] , provide the examples. 
The general form of the four-quark expectation value can be presented as 
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Q12 = (MqTiq • 9 T 2<?K) (i) 

with the subscripts of the matrices including spin, colour and flavour. By using the reduction formula obtained by 
Lehmann, Symanzik and Zimmerman (LSZ) flIT| , p"2"[ we present the matrix elements of Eq.(l) through the expectation 
values in vacuum 

Qn = ^£<0|gT<g.«T i g|0) (2) 

with /„. ss 93 MeV being the pion decay constant. Further calculations are carried out in the framework of a 
factorization hypothesis for the vacuum expectation values. This hypothesis was first formulated by Shifman et al. 
p3| and was advocated recently in ]14| . In this approximation the expectation values of the light quark operators are 
expressed through the well known value [fl5f of the condensate 

(0|qq|0> = (3) 

m u + m d 

with m w<u ,d standing for the masses of the pion and of the light quarks. For example, if "g" is "u" or "d" quark 

w<*rw - --^f^- 

Jit 
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for any state of the pion isotope triplet. 

The negative value of the condensate corresponds to a domination of "disconnected terms" when one of two qq 
pairs comes from the vacuum. The qq pairs in the pion are found to be correlated strongly, with the probability of 
the coexistence of two pairs being much smaller than in the independent pairs picture. 

The four-quarks interactions are presented in Nambu and Jona-Lasinio model [jig] , emerging in its simplest SU(2) 
version. The instanton induced t'Hooft interaction |l7| is often included in the applications nowadays-see e.g. pq . 
Being written in terms of the current quarks, these interactions Vi,2 correspond to the picture "before" the symmetry 
was spontaneously broken. Thus the result of averaging over the pion states (7r|Vi|7r) = (7r| V2 |7r) = is a reasonable 
one. This means that introducing the interactions VJ.,2 does not change the energy of the ground state. Hence, the 
latter remains degenerate until the spontaneous symmetry breaking occurs. 



II. GENERAL EQUATIONS 

We define the expectation value of any operator A in a hadron h as 

(h\A\h) = (h\ f d 3 x \a(x)- (0\A(x)\0)] \h) (5) 



with the vector of state (h\ normalized as (h(k) | h(p)} — 2poS(k — p). Of course, the expectation value (0|A(x)|0) 
does not depend on x. Thus the matrix element {h\ A\h) is defined as the excess of the density of the operator A over 
the vacuum value, integrated over the volume of the hadron. 

Due to the partial conservation of axial current (PCAC) the pion state vector can be expressed through that of the 
vacuum (see, e.g., Q) 

KOr)) = 1 - d,J« a (x)\Q) (6) 

with a being the isospin index while J^ s is the axial current of the light quarks with the corresponding quantum 
numbers: 

J ^( x ) = ^2d c (x) llllb u c {x); J+(x) = ^u c {x) llllb d c {x) ; 

C C 

J °5(X) = Y\ ^(^H'^ 5 ^) Z d c {x)-f^d c {x) ^ 

c ^ 

with "c" being the colour index. 

By applying Eq.(5) for both (tt| and \n) states in the matrix element (7r|^4|7r) of any operator A one finds the 
reduction LSZ formula 

1 

72 

J 7T 

with 

B a = ^ J d 3 xdy dz 5(x - y )5(z - x ) [Q%(z ), [Q%(y ), A(x)}] . (9) 

Here V is the normalization volume, the commutator [X, Y] = XY — YX, while 

QUvo) = J d 3 yJ£ 5 (y) (10) 

is the axial charge, corresponding to the current J" 5 . 

Note that application of Eq.(9) to the quark scalar operator A — uu + dd yields H 

{tt\uu + ddln) - -I— (11) 

m u + md 

with the proper behavior in the chiral limit. It was shown in Jl9[ that the dimensionless value, identified with the 
total number of quarks and antiquarks 

_ (ti\uu + dd\iz) 

can be obtained by using the pion vector of state |fr) with normalization (jr(k)\jr(p)) = <5(k — p). 



(n°\A\n a ) = - 2 (0|B«|0) (8) 
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III. THE FOUR-QUARK OPERATORS 



The general form of the operator, containing the four quark fields q of a specific flavour is A = q^q^q^q^ with 
a, (3, 7, 8 and a, 6, c, d being Lorentz and colour indices. In the applications the product of the fields is time-ordered. 
Thus, considering the product of the field operators at the same space-time point we treat it as the limit of a time- 
ordered product. 

Each of the products of the two-quark fields can be presented as 

C$ = E ^ X 1^ ab - 1 £ qT x yqT x a( X b (13) 

X X,p 

with A p (p — 1 . . . 8) standing for the standard Gell-Mann matrices normalized by the relation Sp A P A T = 2S pT . The 
16 basic Dirac 4x4 matrices are 

r 1 = i , rl = 7m , r 3 = 75 , r£ = 7ai75 , rf = ^{ llllv - 7i/7m ) (14) 

with fj, = 1...4, k = 1...6. 

Thus it is sufficient to study the condensates of the form 

Sxy = (fc|XyrV-? 6 rV|fc) (is) 

a , 6 

and 

R XY = {h\ ]T fr^^/.g^A^/lft) (16) 

a,a' ,b,b' ,p 

with the matrices T(A) acting on Lorentz (colour) indices. 

For the expectation values of the operators, which are antisymmetric in the colour variables 

U XY = (h\ J2 9°rV • q b r Y q b '\h) (S aa , ■ 5 W - S ab , ■ S ba ,) (17) 

a, a' ,b : b' 

we find by using the properties of Gell-Mann matrices 
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Uxy = - Sxy — - Rxy (18) 

with Sxy and Rxy being defined by Eqs. (15), (16). 

The condensates determined by Eqs. (15), (16) contain the quarks of the same flavour. For the case of different 
flavours of quarks: u, d(u, d, s) one can present similar covariant expressions replacing T X,Y by T i ' = T X,Y ■ t l with 
t l acting on the isospin indices. 



IV. EXPECTATION VALUES OF THE LIGHT QUARK OPERATORS 

Here we calculate the expectation values of the 4-quark operators which include "it" and "d" quarks only. The 
commutators which enter Eq.(8) can be calculated by using the relation Jl6| 

{q a a (x),q b p {y)} = ( 7 o)a^ afc( 5(x - y) (19) 

(with {X, Y} = XY + YX) at x = yo, while the anticommutators between all the other quark operators turn to 
zero. Evaluation of the right-hand side (rhs) of Eq.(9) thus leads to the form 

(0\B a \0) = i <0| I d 3 xF a (x)\0) (20) 

with F a {x) being the product of four quark operators. Since the vacuum is uniform, the density F a does not depend 
on x. Hence, we can put 
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B a = F a (0) . (21) 

To illustrate how it works we present the calculation of the operator B, corresponding to the expectation value 
(ir~\uTu ■ uTu\tt~) with T = T x = T Y . The internal commutator in the rhs of Eq.(9) is 

K = J d3 vJ2 [d a (yhol 5 u a (y) , u(x)Tu(x)u(x)Tu(x)] . (22) 



Presenting explicitly d a 7o75W a = J2 7 a ^"(ToTs^c^S an d uTu = J2 @ $ &w^r psu b 5 we employ the fact that the anti 

b 

"0 



commutator of the fields u a a and u b g is given by Eq.(19) while all the others ones are zero. Using also the relation 



{7o, 75} = we obtain 

k = - {d(x)j 5 Tu(x) , u(x)Tu(x)}. (23) 

After similar calculation of the commutator [Q5 , k] we find 

B~ {(uTu) 2 ) = -uTu (uTu + dj 5 Tj 5 d) - dj 5 TuuTj 5 d . (24) 

Here we introduced notation B a (A) for the operator B related to operator A by Eq.(9). 
One can obtain in the same way 

B + ((uTu) 2 ) =B-((uTu) 2 ) (25) 

and 

B° {(uTu) 2 ) = -uTu(uTu + u'y 5 rj 5 u)-^u{'y 5 ,r}u-u{'y 5 ,r}u . (26) 
For the mixed operator A = uTudTd we find 

B ± (uTudTd) = —uTudTd — — {uTuu^T^u + uj^Tddj^Tu + (u <-> d)) , (27) 

while 

B" {uTudTd) = -uTudTd - i (uTudj 5 Tj 5 d + dTdu^T^ 5 u - u{ 75 , T}u d{j 5 , T}d) . (28) 

Of course, the isotope invariance provides the relations B + {(dd) 2 ) = B~ ((flu) 2 ), etc. 
For the important special case T = I Eqs. (24)-(28) take the form 

B^ {(uu) 2 ) = —(uu) 2 — uudd — d"/ 5 uu , y 5 d (29) 

B°((uu) 2 ) = -2 {(uu) 2 + (u l5 u) 2 ) (30) 

and 

B ± {uudd) = -i {{uu + dd) 2 + 2u-y 5 ddj 5 u) , (31) 

while 

B°(uudd) = — 2 {uudd — u/y^udrf^d) . (32) 

Note that Eqs. (24)-(32) are obtained for the operators Sxx defined by Eq.(15), i.e. for the colourless diquarks. 
The operators B(Rxx) with Rxx defined by Eq.(16) are expressed by the same Eqs. (25)-(28) with the matrices T 
being replaced by T p = T\ p with further summation over p. 

By using Eqs. (9) and (18) one can obtain the operator B{A) for the four-quark operator of the general form 

A = ^*f/Pi*-*f^P 2 * (33) 
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with ^ = ( j) being the quark spinor. The operators Pi. 2 are the projection operators in the isospin space. Each of 
them is equal either to P+ or to P_ with P± = (1 ± t^)/2. We defined 

f£ = r fc A" (34) 

for any 4x4 matrix acting on Lorentz indices. 

The matrix A p (p — ... 8) is either the unit matrix (p = 0) or one of the Gell-Mann colour matrices (p = 1, ... 8). 
For the operator A defined by Eq.(33) we obtain 

B a = -i* (75l> a P 2 + T p i15 P 2 t s ) * 
x* ( 7 5rfr Q Pi + rf 75 Pir Q ) * - i ^T P } P 2 ^ 

x (§irfr a r tt Pi* + §rfPir Q r a *) + (1 -» 2) (35) 

with a = +, — , standing for the pion isospin indices. For the specific cases Eq.(35) takes a much more simple form 
since some of the terms vanish: r+P+^P = t_P_\I/ = 0; P_r+\f r = P + T-^f = 0. 

V. FOUR-QUARK OPERATORS CONTAINING HEAVIER QUARKS 

Here we calculate the expectation values of the operators containing heavier quarks. Calculations, similar to those 
described in the previous Section, provide 

B ± (uT x u^T Y ^) = -1 (uT x u + d l5 T x l5 d) ^T Y ^ , (36) 

while 

B° (qT X q^T Y ^) = ~ (qT X q + q^T X l5 q) ^T Y ^ . (37) 

For the operator A consisting of the heavy quarks only, e.g. A = %piT x ipi-4>jT Y 'tpj , we find immediately 

B a {A) = . (38) 



VI. FACTORIZATION APPROXIMATION 

Now we calculate the rhs of Eqs. (24)-(37) in the factorization approximation JD3]. This means that the sum over 
the intermediate states is assumed to be dominated by the vacuum state, i.e. 

(°\Qa,a^,b^,c<lid\ ) = -^(5 l3 5 a p6ab5 k e6 7 s5cd - 
SuS a sS ad 6 jk S^S bc ) (0|qV|0) (0|qV|0) (39) 
with the upper indices standing for the flavour. For the quarks of the same flavour this provides for any 4x4 matrices 

(0\qfPqqfPq\Q) = ^ (s P f ' Sp f> - Sp f p f p ) ((0\qq\0)) 2 (40) 
with T p k defined by Eq.(34). For p = 0, when A p is the unit matrix and T° = Tfc, Eq.(40) takes the form 

{0\qT r qqT s q\0) = ^(Sp T r Sp T s - ^Sp T r T s ) (0\qq\0)f , (41) 

while 
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(0|^qT^T^|0) = -± S P r r r s ((0|w|0)) 2 . (42) 

p=l y 

For the quarks with different flavours i ^ j we find 

(o|*r r ^r s ^|o) = ^ s P r r Spr a <o|ft ft |o)<o|^|o) (43) 



and 



while 



(0|g i r r ^-r a g j |0> = -g^g Spr r r a <0|ft<7i|0)<0|^|0) , (44) 



(o|^r^r^|o) = -\ s P r r r s (o|q l% |o) (o\q jgj \o) , (45) 

which is true for i = j as well, and 



(o|^t^t%|o) - 0. (46) 
P =i 

Assuming also isospin invariance, i.e. (0|uu|0) = (0|Jd|0) = (0|qq|0) we find 

(0|S+ ((uTuf) |0) = ~l (( Sp T) 2 - 1 Sp T 2 ) (<0|^|0>) 2 (47) 
(0\B~ ((OTu) 2 ) |0) - (0|B+ ((ur M ) 2 ) |0) (48) 

while 



(o|s° ((sr«) 2 ) |o) = -1 (( s P r) 2 - 1 s P r 2 - ± s P r 75 r 75 + (s P 75 rf 



.((OI^IO)) 2 (49) 



For the mixed operator A = uTudTd we obtain 



and 



(01^10) = - 1 1 (Sp r) 2 - ^s P r 75 r 75 ) «o|w|o>r (, r >0) 



(o|B°|o) = - I ((s P r) 2 - (s P 75 r) 2 ) . ((o|^|o)) 2 (51) 



We present also the result for the operators containing the coloured diquarks. For the operator J2 P ( u ^^ Pu ) 2 they 
are 

<o|B+|o> = I Sp r 2 ((0|^|0» 2 (52) 
(0|B"|0)= § Sp r 2 «0|^|0)) 2 (53) 

(o|B°|o) - I ( s P r 2 + sp r 75 r 75 ) ((o\q q \o)) 2 , (54) 

while for A = Y Jp uY\PudT\Pd 

(01^10) = ^ s P r 75 r 75 ((o|w|o)) 2 ; (o|s°|o) = . (55) 

Equations (47)-(55) are true for any 4x4 matrix T. If T = T x with T x , being one of the basic matrices defined 
by Eq.(14), the unit matrix T = T 1 is the only one for which Sp T ^ 0. Thus the scalar condensate appears to be 
numerically larger than the other condensates of the form A = (q~T x q) 2 . The values are 
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(0\B ± ((qq) 2 )\0) = -_«0|gg|0)) 2 (56) 
2D 

<0|S°(Cgg) 2 ) |0) = -_((0|^|0» 2 , (57) 



while for A — uudd 
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(OI^IO) = -- ((0|gg|0)) 2 (58) 

(0|i? o |0)= -2«0fe|0)) 2 - (59) 

Deviations of the coefficients in the rhs of Eqs (56)- (58) manifest themselves in the units of the characteristic factor 
1/12-see Eq (41). Also, for the operator A = £ p (uA p u) 2 

(0\B±\0) = - ((OmO)) 2 ; (0\B°\0) = - ((0\qq\0)) 2 , (60) 

while for A = J2 P u\ p ud\Pd 

(01^10)= | «0|gg|Q» 2 ; (0|£°|0) = . (61) 
For the expectation values of the operators containing heavier quarks we obtain for the operator A = qT x q%l)iT Y tpi 
(0|B|0) = ™ (Sp Y x Sp T Y + Sp T x Sp 75^75) (0|g#) (0|^|0> , (62) 
while (0|S|0) = for the operators containing the coloured diquarks. 

VII. SUMMARY 

We solved the problem of expressing the expectation values of four-quark operators in pions through those in 
vacuum. The further specific calculations were carried out in the framework of the factorization hypothesis. 

For the scalar operators of the light quarks the expectation values turn out to be negative. This can be understood 
in the following way. Note that the four-quark condensates take the form of the superposition of the products of the 
two-quark operators f?i,2, i-e., A = 771772 ■ For 771 = 772 = i] Eq.(5) can be presented as 

(h\ v 2 \h) = 2<0|77|0) (h\ V \h) + (h\ ( d 3 x( V (x) - (0|77|0))» , (63) 



if the factorization approximation is assumed. 

The first term in the rhs of Eq.(63) describes the "disconnected contribution" with one of qq pairs coming from 
vacuum. For the expectation values of the operators (qq) 2 , i.e. 77 = qq, the first term in the rhs of Eq.(63) is exactly 
the rhs of Eq.(4). Thus the expectation values of operators (qq) 2 are dominated by "disconnected terms". The same 
refers to the operator A = uudd. For the operators 77 of the form qT x q with X ^ 1 (see Eq.(14)) the "disconnected 
terms " vanish since the scalar condensate is the only one with a nonzero vacuum expectation value. 

For the scalar operators the second term of rhs of Eq.(63) describes the correlation between qq pairs inside the 
pion. It corresponds to deviations of the rhs of Eqs. (56)-(59) from that of Eq.(4). If these pairs would have been 
independent, we would expect the estimation 

(0|B|0) - (-2) • « H?gk ) • f (64) 

to be true. Here n%. being the number of qq pairs of the certain flavour (see Eq.(12)), while V v is the pion volume. 
The latter can be obtained since the pion radius is known to be |2(j] 
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However, the value obtained from Eqs. (56)-(59) appears to be much smaller than the rhs of Eq(64). This means 
that the pairs of light quarks are correlated strongly and that the probability of finding two pairs at the same point 
is much smaller than in the approximation of independent pairs. 

As one can see from the formulae of Sec. 5, for the condensates of the same flavour qT x qqT x q those with X = 1 
are about 10 times larger than the other values of X. In other words, the value of the scalar condensate is about 10 
times larger than the values of the vector, pseudoscalar, pseudovector and tensor condensates. This is because T 1 is 
the only matrix with SpT x ^ 0. Since for each of the matrices T x defined by Eq.(14) |Sp (r x ) 2 | = 4, the condensates 
with 1^1 contain the numerically small factor ^ = For the condensates with different flavours uT x udT x d the 
scalar condensate dominates for the charged pions. For the neutral pions the pseudoscalar condensate is as large as 
the scalar condensate. 

We obtain also the values for the condensates containing heavier quarks. Assuming (0|ss|0) ~ (0|gq|0) with q 
standing for the "it" or "c?" quark, we find the expectation values (Tr\qqss\ir) to be only about twice smaller than 
the value (7r|(gg) 2 |7r). The expectation values, containing two pairs of heavier quarks (s,c, etc.) become zero. The 
expectation values of the four-quark operators in pions in the terms describing the four-quark interactions in NJL 
model are zero. This is the expected result corresponding to the degeneracy of the ground state before the spontaneous 
symmetry breaking. 

We thank M. Eides and T. Gutsche for useful discussions and E.Gcrjuoy for reading the manuscript. Two of us (E. 
G. D.) and (V. A. S.) are grateful for the hospitality of University of Tucbingcn during their visits. The work was 
supported by DFG grant No.438 RUS 113/595/0-1 and by RFBR grant No.00-02-16853. 
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